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Abstract 

The effect of strain on the Landau levels (LLs) spectra in graphene is studied, using 
an effective Dirac-like Hamiltonian which includes the distortion in the Dirac cones, 
anisotropy and spatial-dependence of the Fermi velocity induced by the lattice change 
through a renormalized linear momentum. We propose a geometrical approach to 
obtain the electron’s wave-function and the LLs in graphene from the Sturm-Liouville 
theory, using the minimal substitution method. The coefficients of the renormalized 
linear momentum are htted to the energy bands, which are obtained from a Density 
Functional Theory (DFT) calculation. In particular, we evaluate the case of Dirac 
cones with an ellipsoidal transversal section resulting from uniaxially strained graphene 
along the armchair (AG) and zig-zag (ZZ) directions. We found that uniaxial strain in 
graphene induces a contraction of the LLs spectra for both strain directions. Also, is 
evaluated the contribution of the tilting of Dirac cone axis resulting from the uniaxial 
deformations to the contraction of the LLs spectra. 
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1. Introduction 

Graphene is a two-dimensional material conformed by hexagonal rings of carbon 
atoms with unique physical properties, which make it a material of great scientihc and 
technological interest. Graphene shows an anomalous Hall effect at room temperature. 
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resulting from a linear dependence between the density of electric charge carriers and 
the voltage An extraordinary feature of the electronic properties in graphene 

is the linear relation between energy and momentum of 2pz electrons at low energies 
(less than 1 eV), being different than the usual quadratic energy momentum relations in 
ordinary materials. That dispersion relation [E oc k) causes that the electron transport 
in graphene to be governed by a Dirac-like Hamiltonian, and electrons to behave like 
massless Dirac fermions W- The peculiar electronic properties of graphene have 
consequences on the Landau Levels (LLs) spectra, which are dehned as quantized energy 
states of charged particles in motion under a uniform applied magnetic held, B [6]. 
In an ordinary conductor, the energy of the LLs have a linear dependence with the 
quantization integer n, as {n + 1/2)5. However, in graphene the LLs spectrum is 
quantized according to \/nB [THU]. 

In recent years, many works have been devoted to hnd practical methods to control 
graphene’s properties. A possible way to modulate the electronic, vibrational and 
transport properties is by performing a deformation on the graphene sample or 

bilayer graphene m- Thus, actual studies identify three main effects of uniaxial strain 
on low-energy band structure in graphene [131 - 120] . First, a slipping of the Dirac points 
out of the high symmetry points K and K’. Second, a distortion in the traversal section 
of the Dirac cones, breaking the isotropy of the Fermi velocity; and third, a small 
vertical axis tilting of Dirac cones, which could be negligible in most of the situations. 
Considering that the LLs in graphene only depend on the shape of Dirac cone cross 
section, the second effect is expected to have more influence on the structure of LLs 
spectra. As we show here, the tilting of Dirac cone axis also has a non negligiable effect 
on LLs for high strains (15-20%). 

Some theoretical works use an anisotropic mass model to hnd the LLs in graphene 
through topological considerations HI. Another authors frequently analyze the prop¬ 
erties of strained graphene using the Tight-Binding (TB) approximation, including the 
effect of deformation on the atomic distances through the scaling of the hopping param¬ 
eters Hal HMHl [2IH29]. This hopping renormalization is commonly modeled with an 
exponential decay |13] or using the Harrison’s scaling rule Vpp-jr (x 1/P [Hj. Both renor- 
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malizations could fail beyond the linear elastic regime since the Poisson ratio changes 
with the strain. Indeed they could have a different dependence for each strain-type 
considered. An accurate and precise way to obtain the hopping parameters in strained 
graphene is through a htting of the TB Hamiltonian to the energy bands obtained from 
Density Functional Theory (DFT) calculations [30] or from experimental data. 

Instead of a hopping renormalization, an alternative way is proposed in the present 
work. We use a renormalized linear momentum in an effective Dirac-like Hamiltonian. 
The coefficients of the linear momentum are related with the geometrical parameters 
of the distorted Dirac cone, and can be calculated from a htting to the energy bands 
obtained with a DFT calculation. Then, we apply a minimal substitution in the free- 
held ehective Dirac-like Hamiltonian to get the LLs spectra. In particular, for uniaxial 
strain we found that the LLs spectra is contracted as a function of the deformation 
along the Zig-Zag (ZZ) and Arm-Chair (AC) directions. This contraction of the LLs 
spectra is due to the renormalization of the Fermi velocity with the strain, which is 
reduced by the stretching along these directions. In addition, we have evaluated the 
contribution of the tilting of Dirac cone axis to the contraction of the LLs in uniaxially 
deformed graphene 

2. Landau levels in graphene 

To clarify our methodology, we show how to obtain the LLs in unstrained pristine 
graphene from an effective Dirac-like Hamiltonian in a low energy and magnetic held 
regime. We treat the dynamics of an electron moving in a graphene sheet under a 
uniform magnetic held B = Bz perpendicular to the direction of propagation. We 
consider the electrons in graphene as massless Dirac fermions, having a linear dispersion 
relation (£' oc fc) [3]. With H = 0, the Hamiltonian H only depends on the linear 
momentum p and can be represented by a 4 x 4 block diagonal matrix [ISlEs] 



( 1 ) 


where vp is the Fermi velocity and a are the Pauli matrices acting on the pseudo-spin 
space, which discriminate between the contribution of the two triangular sublattices 
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present in graphene. Each block in the Hamiltonian ([^ represents the K and K' 
valleys, coinciding with the high symmetry points in absence of strain. For H 7 ^ 0, 
we do the minimal substitution p —)■ p + eH in the free-held Hamiltonian Q with the 
Landau gauge A = xBy pin]. Due to the Py conservation, the electron’s wave-function 
in graphene can be expressed using variables separation, having the form 

^ (x, p) = (2) 

where ^(x,p) and v{x) = {(p^~^\x),(l)^~\x)) are four-component vector functions with 
= {f^\x),g^j^\x)) and \x),f^ describing the pseudospin 

({/^(x), pB(a^)}) and the valley {K or K’) states ({±}). 

We focus in the K valley, substituting (|^ and ([^ in the Dirac equation iL^(x, y) = 
E\[^(x, p), we obtain 

Vf {Pox - iPoyiky)} ^ 2 ^ 

Vf {Pox + iPoy{ky)} fA^\x) = Eg^^\x), 

with Pox = Px and Poy{ky) = hky -|- exB. Decoupling the system of equations ([^ and 
using the commutator \px)Poy{ky)] = —iheB, we get 

vHpI + {hky + exH)2}/W(x) = {E^ - vlheB}f^/\x), (4) 

being similar to the quantum harmonic oscillator equation. Therefore, for Dirac elec¬ 
trons in the presence of a uniform perpendicular magnetic held, the LLs spectra is given 

by [inilll] 

En = sgn{n)hw^^/\n\, = vf^‘^^, ( 5 ) 

with n = 0, ±1, ± 2 , ...and degeneracy of 4:SB/{h/e), where S is the sample area 
and {h/e) is the magnetic hux quantum. The four-fold degeneracy of LLs in deformed 
graphene (n 7 ^ 0 ) is due to a two-fold pseudospin and a two-fold valley degeneracy. 
While that for n = 0, there is a two-fold degeneracy because the valley index is the 
same as the sublattice index. This result is diherent from the LLs spectra for con¬ 
ventional conductors = hu{n + |), where each level has a constant separation, 
while in graphene the LLs separation (|^ decreases as |n| increases. This behaviour has 
been conhrmed experimentally mu [To] . It is important to mention that htting the 
expression ([^ to the experimental spectrum, the value oi vf can be obtained da- 
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3. Landau levels in strained graphene 


It is known that the Dirac cones lose their isotropy when graphene is under non¬ 
isotropic strain and consequently, the effective Dirac Hamiltonian of Q is not valid 
anymore. To overcome this problem and thinking in a general case, i.e. an inhomoge¬ 
neous strain [22H2U 1311 - 135] , we propose a renomalization in the linear momentum p for 
the Hamiltonian ([^ 

p = a{x, y)p^x + 6(x, y)pyy, (6) 

where a(x, y) and h{x, y) are dimensionless functions. In pristine graphene, the Dirac 
fermions move with an approximate vp of 10® m/s [2]. For strained graphene, the Fermi 
velocity is anisotropic and space-dependent [221123]. The present approach implicitly 
includes the anisotropy and spatial dependence of the Fermi velocity through the di¬ 
mensionless functions a{x,y) and b{x,y). Substituting 0 in the Hamiltonian Q for 
the iF valley 

/ 0 p^a{x,y) -ipyb{x,y) 

H = vf \ 

\ a{x,y)p^ + ib{x,y)py 0 

with a vector function T(x, ?/) = describing the pseudospin {{fA,gB})- 

For a nonuniform magnetic or pseudomagnetic field perpendicular to the strained 
graphene sheet |32HS1|, we do the minimal substitution p ^ = p + eA^ in the 

free-held Hamiltonian Q, where A{x,y) is a potential vector that can depend of the 
inhomogeneous strain and acts different in to each valley [3211M] . We focus in the K 
valley, and decoupling the 2x2 equation system obtained from the Dirac equation 
= ET, we get 

{(a7r+)2 + (67r+)2 - Ai(a7r+ + i&7r+) + i[a7r+, 67r+]|/j+^(a:, ?/) = ^ff\x,y) 

F 

|(a7r+)2 + (&7r+)2 - Ai(a7r+ + ibn^) - i[a7r+, &7r+] + X2jg^B\x,y) = -^g^p\x,y) (8) 

where (a7r+)^ = {apxY + 2ea?‘A'^‘^px + — ieha-^{aA^), a similar expression is 

obtained for doing the changes x ^ y and a ^ b. The commutator is given 

by [a7r+,67r+] = -ih{a^){py + eA':!^) + ih{b^){px + eA'^) - iheabB^ with B^{x,y) the 
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magnetic field generated by A~^, the Xi{x,y) = and X 2 {x,y) = h?{—a^ — 

— b-^^)}. For the K' valley, the same decoupled equation system is 
obtained with g^^\x^y) and f^\x,y) instead of f^\x,y) and < 7 ^^(a:,?/) respectively, 
with superscripts (—) in the tt operators and A vector potential. The 2x2 partial 
differential equation system (|^ is an alternative way to describe the electron dynamics 
in strained graphene with deformed Dirac cones and in the presence of a nonuniform 
magnetic held from the anisotropy parameters with a spatial dependence. If we know 
the exact form of a{x, y) and h{x, y) from the dispersion relation or the Fermi velocity, it 
is possible to obtain the LLs solving the equation system (|^. When an inhomogeneous 
unidirectional strain is applied, these functions must depend only on x (y) due to the 
translational symmetry on y (x). Assuming x dependence for the anisotropy parameters 
and using the Landau gauge A = xBy, the Py component of the linear momentum is 
conserved. Thus, the equation system in (|^ is reduced to 



where R{x) = 2b{x) ^ equations (|^ are a general differ¬ 

ential equations for the electron’s wave function in nonuniform unidirectional strained 
graphene, which can be solved using the Sturm-Liouville theory |1^ with eigenvalues 
Xn = El/h^v%. 

In previous studies on uniaxially strained graphene [IMS ED 122], it was shown that 
the Dirac cones have an elliptical cross section for low energy regime |i?| < 1 eV. As we 
shown in the Figure the semi-major axis of the elliptical cross section is always along 
the tension direction. Hence, the functions a{x) and h{x) are approximately constants 
and they are related to the geometrical parameters of the ellipse, such as the semi¬ 
major axis (A = E/ahvp for ZZ and A = E/bhvp for AC direction) and semi-minor 
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axis {C = E/bhvp for ZZ and C = E/ahvp for AC direction), or equivalently, scaling 
the energy axis with q = E/hvp, then a = cot a = q/k^ and b = cot (3 = q/ky, where 
a and {3 are the extremal elliptical cone angles. For pristine graphene a = (3 = tt/A. 
With a and b independent of x, the equation ([^ is reduced to the quantum harmonic 
oscillator equation. In the present case we obtain that the LLs spectra is given by 


En = sgn{n)\^hw^ \/\n\, = VF\j —— 


( 10 ) 


with n = 0, ±1, ±2, ... and degeneracy of AS'B/{h/e), where S' is the deformed sample 


area. From equation (10), we can see that the strained graphene LLs depend on the 
geometrical parameters of the deformed Dirac cone, the respective cyclotron frequency 
and the quantum number n. Since the expressions ([^ and (10) are similar, we 
have dehned ^ as a parameter to measure the contraction < 1) or expansion 

> 1) of the LLs spectra under the same magnetic field. It should be noted that 
in others works, the concept of renormalized Fermi velocity v*f is commonly used na, 
and which is related with our parameter hy v*p = It is important to mention that 
when the cross section of the Dirac cone does not have an exact elliptical form, it is 
necessary to solve the differential equation (|^ to obtain the spectrum of LLs. In our 
best knowledge, there are not experimental reports on the effect of strain on the LLs 
spectra of graphene. 


4. Computational details 

In order to know the geometrical parameters of strained graphene Dirac cones, DFT 
calculations were performed within the plane-wave pseudopotential framework, as im¬ 
plemented in the Quantum ESPRESSO (QE) package |16]. The exchange-correlation 
functional was treated with the PBE parametrization of the generalized gradient ap¬ 
proximation m- Core electrons were replaced by an ultrasoft pseudopotential from the 
QE-PSLIB database [IB], and valence wave functions (charge density) were expanded 
in plane waves with a kinetic energy cut-off of 40 Ry (320 Ry). 

We considered uniaxial strain up to 20% of deformation along the AC and ZZ 
directions (see onsets in Figure [^. Because such amount of strain is beyond of the 
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Figure 1: Contour plots of the Dirac cones in pristine (center) and uniaxial strained (left and right) 
graphene. The uniaxial strain is along the ZZ (left) and AC (right) directions for a deformation of 
10% for both cases. The origin of both axis corresponds to the Dirac point; note that the strain shifts 
the Dirac point from the K point. The cross marks refer the positions of the high symmetry K point 
in deformed graphene. The numbers on each curve indicates the value of the corresponding energy in 
eV. 


elastic regime of graphene, for a given strain, the Poisson ratio was obtained by a direct 
minimization of the electronic total energy. In each one of the steps to compnte the 
Poisson ratio, we relax the carbons positions nntil the interatomic forces were 0.0001 
Ry/Bohr or less. 

For structural and energetic calculations, we employ a grid of 36x36x1 /c-points 
within a 0.01 Ry of cold smearing [38]. The E{kx, ky) surface was interpolated from a 
denser grid of 144x144x1 fc-points and 0.001 Ry of smearing. In all cases, we left 10 
A of vacuum space between successive layers to avoid spurious supercell interactions. 

As expected, the Dirac cone cross section has an approximately elliptical shape, 
as shown Figure We notice that the slipping of the Dirac points out of the high 
symmetry points K and K’ becomes evident for a deformation of 10%, for both de¬ 
formation directions. In order to obtain the geometrical parameters of the deformed 
Dirac cones, we have htted the conduction band around a Dirac point with an energy- 
cutoff of 0.3 eV for a range of deformation up to 20%, using the dispersion relation 
E = hvp-sja^k‘1 -b b'^ky, where a and b are the fitting parameters and {k^, ky) the linear 
momentum coordinates around the Dirac point. In our elliptical cone approximation, 
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Figure 2: (a) Geometrical parameters a and b of the Dirac cone as a function of tensile strain 
(e) along the AC and ZZ directions, (b) Evolution of the ^ parameter for strained graphene as 
a function of tensile strain (e) along the AC (red line) and ZZ (blue line) directions with and 
without tilting of the Dirac cone axis. The black curve correspond to the renormalized Fermi 
velocity obtained from a TB approximation considering interactions up to second nearest 
neighbors in deformed graphene by Goerbig et al. m- 


the values for the standard deviation are less than 0.003 eV, proving the accuracy of 
this approximation. 


5. Results and Discussion 

We found that for both strain directions (AC and ZZ), in the whole range of studied 
deformations, the LLs spectra is contracted (.^ < 1) with respect to the pristine case 
(^ = 1) under the same magnetic field. It is interesting to note that the value of 
and hence the distance between LLs, decrease as the uniaxial deformation is increased. 
This behaviour can be explained in terms of the cyclotron orbit motion change induced 
by strain. Thus, if graphene is stretched, the cyclotron orbit motion has a mean-radius 
bigger than the equilibrium case, then the value of total energy decreases. From this 
perspective, the expansion case in the LLs spectra may occur when graphene sample is 
contracted. 
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In Figure]^ we show the plots of a, b and ^ for different values of deformation (e) 
along the AC and ZZ directions. From the Figure |^a) we observe that for e < 10%, 
a {b) in ZZ (AC) has approximately the same values than b (a) in AC (ZZ). Thus, 
for a hxed value of deformation, the htted cones for AC and ZZ are practically the 
same but one rotated with respect to the other by 90°. Therefore, the effect of AC 
and ZZ uniaxial strains on LLs is expected to be practically the same for strain up to 
10% as we can see in the Figure |^b). This point can be explained if we consider that 
the parameter ^ represents the ratio between the renormalized Fermi velocity and the 
Fermi velocity in pristine graphene. The renormalized Fermi velocity can be seen as an 
effective velocity of the anisotropic Fermi velocity in the cyclotron motion. Considering 
that the cones have identical shape for AC and ZZ directions up to 10%, we expect 
similar values of the renormalized Fermi velocity and similar contractions in the LLs 
spectra for both directions. Also, we observe that ^ has a nearly linear behaviour in 
this range of deformations, with ^ ^ 1 — 0.7e. This result suggest that the degree of 
deformation (e) in a graphene sample can be estimated by extracting the value of ^ 
from the LLs spectra. 

To evaluate the contribution of the tilting of the Dirac cone axis to the contraction of 
the LLs in uniaxially deformed graphene, we have htted the conduction band using the 
dispersion relation of the form E/hvp = aokx + boky + ■^/a'^kl + b'^ky, which is obtained 
diagonalizing the Weyl Hamiltonian [13] , where Oq and bo are responsible of the tilting 
of the Dirac cone axis. With an appropriate rotation of our system, we can use the 
expression ^ = \/~^{l — al/a? — b‘l/b‘^Y^^ [15], instead of ^ = v^. In all cases we obtain 
bo ~ 0, having a tilting only in the kx-Bxis, in agreement with others authors [15]. The 
tilt can be neglected for strain up to 15% as we can see in the Figure [^b), where the 
calculations of ^ with and without tilting are shown. For strains along the ZZ direction, 
the tilt could have observable effects on the LLs spectra for deformations larger than 
15%. In order to compare the present results for ^ with previous calculations reported 
in the literature, in Figure |^b) we include the result obtained from a TB effective 
model reported in [15], where a Harrison’s scaling rule for the hopping parameters 
was used. We can see that the TB approximation, predicts a descreasing behaviour of 
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the renormalized Fermi velocity as a function of strain, in qualitative agreement with 
the DFT calculations, but there is an important quantitative difference, which can be 
attributed to the relaxation of the carbons positions and the Poisson ratio changes 
beyond of the elastic regime included in the DFT calculations. The contraction of the 


(a) 



Tensile strain e (%) 


(b) 
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Figure 3: Evolution of the LLs spectra as a function of tensile strain along the (a) AC (red) 
and (b) ZZ direction (blue). The label in each curve corresponds to the Landau level (n). 


LLs spectra as a function of the tensile strain is shown in the Figure for both strain 
directions (AC and ZZ). As we can see, the contraction of the LLs spectra for the ZZ 
deformation is larger than for the AC deformation. The LLs spectra does not present 
lifting of the two-fold valley degeneracy since our prediction is based in the low energy 
and low magnetic field regime [T^ |28] . 

In order to show the effect of strain on the LLs Density of State (DOS) under a 
uniform magnetic field, we have calculated the DOS for both strain directions with a 
deformation of 20%, and comparing them with the corresponding to pristine graphene 
under the same magnetic field. The DOS for strained graphene along the AC and 
ZZ directions are shown in the Figures |^a) and|^b), respectively. In these plots, we 
observe that the contraction of the LLs spectrum is larger along the ZZ than for AC 
direction, because ^zz = 0.745 and ^ac = 0.801 for e = 20%. The DOS contraction 
with respect to equilibrium spectrum is a result of the LLs spacing reduction due to the 
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Figure 4: Calculated LLs density of states (DOS) for (a) pristine graphene (gray) and uniaxially 
strained along the AC direction (red) for deformation of 20% {^ac = 0.801) and (b) pristine 
graphene (gray) and uniaxially strained along the ZZ direction (blue) for deformation of 20% 
i^zz = 0.745). 


modulation of the Fermi velocity by the strain. Thus, the present results show clearly 
that the LLs spectra in graphene can be modulated via uniaxial strain. 


6. Summary and final remarks 

We have proposed a geometrical approach to consider the deformation of the Dirac 
cones in strained graphene introducing a renormalized linear momentum in the effec¬ 
tive Dirac Hamiltonian for massless fermions. We found an analytical expression for 
the energy spectrum of the LLs, which is a function of the Dirac cones deformation. 
In particular, we found that uniaxial deformation in graphene induces contraction of 
the LLs spectra for both AC and ZZ strain directions. The present hndings help to 
set the Landau levels spectroscopy in strained graphene, which could be used for mea¬ 
sure quantities of interest as the anisotropic Fermi velocity. Hall resistance and related 
electronic properties [5]. The present model offers a simple way to relate electronic, 
vibrational and transport properties with the geometry of the anisotropic Dirac cones. 
Finally, it is important to note that the present approach can be applied to analyze 
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other two-dimensional materials beyond graphene presenting Dirac cones |1]. 
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